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ABSTRACT 


This  report  presents  the  first  semi-empirical  sol¬ 
ution  suggested  to  describe  the  traction  versus  slip 
relationship  for  rigid  wheels. 

The  kinematics  of  a  slipping  wheel  are  analyzed  and 
the  results  are  used  to  express  the  shear  deformation  as 
a  function  of  slip  which  in  turn  is  utilized  in  a  new  shear 
stress-strain  relationship.  The  shear  stresses  are  express¬ 
ed  at  every  point  of  the  tire-soil  interface  surface  and 
integrated  to  obtain  the  traction. 

Other  established  information  is  also  briefly  dis¬ 
cussed,  in  order  to  complete  the  description  of  the  state- 
of-the-art. 
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DEFINITIONS,  SYMBOLS 


R  Motion  Resistance  (Lbs.) 

N  Normal  Component  of  the  Resultant  Soil  Reaction  (Lbs.) 

T  Tangential  Component  of  the  Resultant  Soil  Reaction  (Lbs.) 
q(c  Central  Angle  Associated  with  the  Resultant  Soil  Reaction  (o) 
W  Vertical  Load  on  the  Wheel  Axis  (Lbs.) 

D  Diameter  of  the  Wheel  (in.) 

Moment  Due  to  Bearing  Friction  (inch  Lbs.) 

£>  Factor  of  Bearing  Friction 

z  Vertical  Distance  between  a  Point  on  the  Wheel  Perimeter 
and  the  Undisturbed  Ground  Level  (in.) 

Sinkage  of  the  Wheel  (in.) 

Central  Angle  Associated  with  Zq  (o) 

p  Pressure  (stress)  Normal  to  the  Boundary  (Lbs. /in. ^) 

k  Cohesive  Modulus  of  Sinkage  (Lbs/in.  n+ 1 ) 

c 

kj  Frictional  Modulus  of  Sinkage  (Lbs/in. n+^) 

n  Exponent  of  Sinkage 

C/C  Central  Angle  Associated  with  Z 

0  2  rr  - 

M  Net  Torque  Applied  to  Drive  the  Wheel  (inch  Lbs.) 

DP  Drawbar  Pull,  Net  Traction  (Lbs.) 

Vs  Velocity  of  Slipping  (in. /sec.) 
io  Slip 

V^  Theoretical  Velocity  of  the  Wheel  (in. /sec) 
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u  Angular  Velocity  of  the  Wheel  (Rad/sec) 

v  Horizontal  Component  of  the  Velocity  Vector  at  A  Point 
x  on  the  Wheel  Perimeter  Associated  with  V  =  ut  (iiu/sec) 

v  Vertical  Component  of  the  Velocity  Vector  at  a  Point  on 
y  the  Wheel  Perimeter  Associated  with  V  =  ut(io/sec) 

C  Instantaneous  Center  of  Motion 

r  Rolling  Radius  (in.) 

y.  Ordinate  of  C  (in.) 

e 

I  Unit  Vector  in  X  Direction  (iru/sec) 

J  Unit  Vector  in  y  Direction  (in/scc) 

P  Angle  Enclosed  by  the  Velocity  Vector  and  the  X  Axis  (o) 

0  Central  Angle  Associated  with  the  Point  of  the  Wheel 
v  at  which  the  Velocity  Vector  is  Vertical. 

-  2 
a  Acceleration  Vector  (iit/sec  ) 

ip  Radius  of  Curvature  (in.) 

o 

s_  Shear  Strength  of  the  Soil  (Lbs. /in.  ) 
max 

2 

c  Cohesion  (Lbs. /in.  ) 

4  Angle  of  Internal  Friction 

Angle  of  Failure  (o) 

2 

s  Shear  Stress  (Lbs, /in.  ) 

j  Shear  Deformation  (in.) 

K  Tangent  Modulus  of  a  Shear  Stress  Strain  Curve  (in.) 
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INTRODUCTION 


The  problems  related  to  rigid  wheel  behavior  in  soft 
soil  were  first  analyzed  by  Bernstein  He  derived  re¬ 

lationships  between  the  load  on  a  towed  wheel  of  given  geom¬ 
etry  and  the  sinkage  as  well  as  the  motion  resistance.  The 
basic  assumption  employed  by  Bernstein  was  that  the  normal 
soil  stress  acting  on  the  wheel  surface  is  proportional  to 
the  square  root  of  the  vertical  distance  of  the  point  in 
question  and  the  ground  level.  Thus,  the  influence  of  the 
soil  on  the  stress  distribution  was  accounted  for  by  a  factor 
of  proportionality  (k).  Russian  agricultural  engineers  have 
introduced  another  '"soil  characteristic"  by  using  a  general 
exponent  (n)  instead  of  the  square  root.  .  Most  of  the 
Russian  publications,  however,  follow  Gerstner  by  assuming 
n  =  1,  to  avoid  mathematical  complications  which  arise  from 
the  non-integer  exponent.  In  1950,  Garbari  ,  published  a 
very  useful  paper  in  which  he  suggested  some  new  notions  such 
as  the  "critical  inflation  pressure". 

A  further  development  in  the  adaptability  of  the  empir¬ 
ical  pressure  distribution  equation  was  accomplished  by 
( 4 ) 

Bekker  ,  who  suggested  a  new  factor  of  proportionality  which 
is  practically  independent  of  the  size  of  the  penetrometer 
footing,  or  that  of  the  wheel-soil  interface  if  the  aspect 


ratio  is  not  allowed  to  be  less  than  5:1.  Bekker  rederived 

(5) 

Bernstein's  equations  using  the  new  soil  value  system 

Vincent  and  more  recently  Hcgedus  observed 

some  new  relationships  as  to  the  normal  pressure  distribu¬ 
tion  under  a  wheel,  but  no  method  has  been  devised  as  yet  to 

describe  it  analytically. 

1 fl 1  (91 

Tanaka  \  and  Phillips  ,  analyzed  the  equilibrium 
of  wheels  in  their  presentations  at  the  First  International 
Conference  on  the  Mechanics  of  Soil-Vehicle  Systems.  Both 
of  the  authors  included  the  tangential  forces,  hitherto 
neglected. 

Uffelmann  conducted  tests  by  means  of  a  special 

apparatus  enabling  him  to  measure  the  tangential  "rim 
stresses"  and  suggested  a  simple  plastic  theory  of  rut  forma¬ 
tion  providing  an  estimate  of  rolling  resistance.  His  equa¬ 
tions  are  equal  to  Bekker' s  equations  for  zero  exponent,  or 
in  other  words,  they  are  valid  for  bearing  capacity  estima¬ 
tions  (zero  sinkage).  Therefore,  no  conclusions  as  to  the 
sinkage  may  be  obtained  on  this  basis. 

An  excellent  analysis  of  the  equilibrium  of  wheels  by 
Schuring  also  stresses  the  importance  of  the  tangential 
forces  In  addition,  it  dwells  on  the  possibility  of 

using  dimensionless  analysis  techniques  for  the  wheel  pro¬ 
blem. 
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The  latter  approach  lias  also  been  elaborated  on  by 
Hicks  and  Nuttal  at  the  First  International  Conference 
(12  and  13). 

The  kinematics  of  a  rigid  wheel  and  that  of  a  soil 
particle  were  analyzed  by  Andreyev^l4\  and  other  Russian 
investigators  in  great  detail.  .  Their  basic  assumption, 
however,  that  the  displacement  of  a  soil  particle  and  the 
force  acting  on  it  are  collinear  is  incorrect  and  leads  to 
certain  ranges  on  the  wheel  perimeter  according  to  slip 
conditions  which  are  questionable  on  the  basis  of  test  re¬ 
sults  obtained  by  Vincent  and  Uffclmanh./  As  for  soft  tires, 

Omelyanov's  empirical  equation  for  resistance  prediction  is 

(15) 

frequently  published  in  Russian  textbooks  , 

Another  Russian  paper  by  Ageykin  reveals  an  equa¬ 

tion  for  sinkago  similar  to  that  published  by  Bekker  and  the 
writer  <17'  l8'  l9'  20). 

(21) 

Experiments  by  the  National  Tillage  Laboratory  ,  by 

the  Waterways  Experiment  Station  of  the  Corps  of  Engineers, 
(22)  (23) 

U.  S.  Army  ,  by  Soehne  ,  and  by  a  group  consisting 

of  researchers  from  the  National  Tillage  Laboratory  and 
Michigan  State  University  ^4)^  elucidated  some  important  as¬ 
pects  of  the  pressure  distribution  under  soft  tires  operating 
in  deformable  soils. 
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Kerr  and  Chapoux  carried  out  extensive 

research  projects  on  tire  behavior  in  sand.  Weinblum  and 
(21) 

Orlowski  ,  as  well  as  other  Israelian  researchers  are 

also  engaged  in  investigating  wheels  in  sandy  soils. 

They  found  that  an  exponential  expression  is  suitable  to 
describe  the  traction  slip  relationship.  Soehne  and  Sonnen 

/nQ\ 

arrived  at  the  same  approximation  in  their  broad  ex¬ 
perimental  studies  presented  in  Turin.  They  emphasized  the 
difficulties  which  arose  with  the  application  of  the  Bekker 
soil  value  system  in  multi-layer  soils. 

Turnbull  and  Freitag  accounted  for  a  massive  tire  test 

program  which  is  being  conducted  at  the  Waterways  Experiment 
(29) 

Station  .  The  authors  investigate  the  change  in  per¬ 

formance  due  to  multiple  passes,  beside  the  traction-slip 
problem. 

A  new,  radial  ply,  tire  design,  introduced  by  the  Pirelli 

Company  ,  proved  to  be  superior  to  the  conventional  one, 

according  to  tests  conducted  by  the  Company  itself,  by  the 

National  Tillage  Laboratory  an(j  by  tbe  por(j  Motor  Co/32). 

Bekker's  condual  tire  concept  appeared  to  be  promising, 

in  view  of  theoretical  considerations  and  some  limited  test 
(3*1) 

results  .  Richey  suggested  the  mounting  of  a  radial  ply 

tire  on  a  narrow  rim.  The  idea  proved  to  be  advantageous 

(32) 

over  a  conventional  rim  equipped  with  a  radial  ply  tire 
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Various  aspects  of  tire  behavior  on  hard  ground  have 
been  analyzed  by  a  great  number  of  researchers  and  in¬ 
stitutions.  An  excellent,  if  not  quite  up-to-date,  digest 
of  these  studies  has  been  published  by  Hadekel*34*. 
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OBJECT 


The  purpose  of  this  research  project  was  the  establish¬ 
ment  of  simple  equations  based  on  empirical  soil  stress- 
strain  relationships  which  would  enable  one  to  predict  the 
sinkage,  the  motion  resistance,  and  the  traction  of  a  rigid 
wheel  and  a  soft  tire.  While  there  have  been  empirical  sol¬ 
utions  suggested  for  the  first  two  problems,  the  question  of 
traction  of  wheels  has  not  been  evolved  farther  than  the 
statement  that  the  tractive  force  is  the  sum  of  the  horizon¬ 
tal  components  of  the  stresses  acting  along  the  soil-wheel 
interface  surface.  A  number  of  researchers  have  stated  that 
a  rigorous  solution  represents  an  impossible  task  at  present. 
This  i 8  probably  true,  since  the  basic  theory  of  soil  mechan¬ 
ics  has  failed  to  answer  such  basic  questions  as  the  true 
stress-strain  relationship  prevailing  under  a  plate.  No 
wonder  that  a  solution  for  the  involved  three-dimensional 
problem  of  a  tire-soil  interaction,  in  which  the  behavior  of 
the  carcass,  and  the  effect  of  the  lugs  represent  even 
further  obstacles  in  the  analysis,  has  not  been  considered 
obtainable.  Nevertheless  it  is  felt  that  adequate  empirical 
knowledge  is  available  on  soil  behavior  to  make  our  goal 
feasible. 


6 


SUMMARY 


Exact  equations  of  equilibrium  for  driven  wheels  oper¬ 
ating  in  soft  soil  are  presented  for  which  a  formula  for 
motion  resistance  is  derived. 

Next  a  detailed  analysis  of  the  kinematics  of  a  slip¬ 
ping  wheel  is  presented  and  a  geometric  representation  of 
the  velocity  distribution  is  introduced. 

An  empirical  soil  shear  stress-deformation  equation 
is  proposed.  This  equation  is  then  combined  with  the  re¬ 
sults  of  the  previous  analysis  and  an  equation  for  the  gross 
traction  of  a  tire  is  derived. 
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CONCLUSIONS 


The  theoretical  investigation  leads  to  the  following 
conclusions! 

The  perimeter  of  the  wheel  may  generally  be  divided 
into  four  regions  depending  on  the  direction  of  the  velocity 
vectors  that  belong  to  the  pertinent  points. 

When  the  sinkage  is  such  that  the  region  characterized 
by  vectors  pointing  forward  and  downward  is  a  part  of  the 
wheel-soil  interface  (in  case  of  driven  wheels)  the  tan¬ 
gential  stresses  result  in  resisting  forces.  Hence,  optimum 
traction  is  reached  when  the  sinkage  is  not  greater  than 
the  height  of  the  point  to  which  a  vertical  velocity  vector 
belongs. 
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RECOMMENDATIONS 


It  is  recommended  that  the  theoretical  method  presented 
in  this  paper  be  thoroughly  checked  by  a  series  of  experi¬ 
ments  to  establish  its  applicability  and  its  limitations 
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THEORETICAL  ANALYSIS 


The  tire  may  be  considered  rigid  if  the  inflation 
pressure  is  higher  than  the  so-called  critical  pressure 
(defined  in  Part  II), 

To  express  the  sinkage  and  the  motion  resistance  of  a 
towed  wheel  operating  in  soft  ground,  consider  Figure  U 


It  is  assumed  that  the  problem  is  two  dimensional. 

The  equations  of  equilibrium  are  the  following: 

R  -  N  sin  o<  +  T  cos  c*.  =0  (1-a) 

0  o 

W  +  T  sin  cXQ  -  N  cos  c<  =  0  (1-b) 

T  |  -  Mfa  =  0  . (1-c) 

where  naturally. 

Thus,  the  magnitude  of  T  can  be  found  if  W  and  9  are 
known.  (Equation  1-c).  The  tangential  component  of  the  re¬ 
sultant  is  a  function  of  slip.  This  relationship  will  be 
derived  on  page  38.  Thus,  it  is  theoretically  possible 

to  establish  the  magnitude  of  the  negative  slip  which  occurs 


Figure  2. 
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Equation  (1)  may  be  rewritten  as  follows: 

dN  sinoC  =  0  ...  (2a) 


i  may  oe 


R  -  d  T  co s/ 

o 


i  a  o  x  u 

./• 


-/ 

i  /• 


o  fro 

dN  cos  dt  +  y'  d  T  sin  dC,- 


0  .  .  .  (2b) 


d  T  -  M.  =  0 

b 


Thus,  if  dT  =  f  (oC)  and  dN  =  g  (®0  is  known  the  sinkage  and 
the  motion  resistance  may  be  found. 

Note  that: 


=  2  U  '  cos^0> 


(3) 


Bekkers  approximative  equation  will  be  used  throughout 
this  paper  for: 

p  =  $[  =  f(z)  =  gK)  =  h  (9) 


p  - 


(\  +  %)J(cos  0  -  cos|Jo)  | ] 

V  . (4) 


Fig.  3 
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If  one  assumes  that  is  negligible,  then  T  =  0 


R  =  dN  sin  = 


2tt  - 


_dN  sin  0  .  .  .  .  (5) 


Using  equation  (4): 


"  =  I2,  =  (IT  +  kf)  "  {§(c°s  8  -  cos  o)J  n  Sin  8 


=  *><r+ k^>  (f)n+1  f1  - cos  P  oj n+i 


From  equation  (4): 


1  -  cos 


Thus: 


k  +  b  k 


n  +  1 


Equation  6  was  first  derived  by  Bekker1  ,  using  energy 
considerations.  It  enables  one  to  approximate  the  rolling 
resistance  of  a  towed  wheel  for  small  slippages  only  because 
the  tangential  forces  have  been  neglected. 
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Bekker  and  Hegedus  have  attempted  to  remedy  this  short¬ 
coming  by  introducing  the  so-called  bulldozing  resistance. 
(5)  (37). 

Equation  (2b)  cannot  be  solved  in  a  closed  form  even  if 

one  assumes  tiT  =  0  for  all  A.  Dekker  has  performed  the  in- 

(5 )  * 

tegration  by  neglecting  the  difference  between  x  and  x  , 

(Fig.  4)  and  by  considering  the  first  two  terms  in  the  bi¬ 
nominal  series  of  n 

[l  -  (Z0  -  z>] 

Ehrlich  (36)  improved  the  accuracy  of  Bekker's  solution 
by  considering  three  terms  in  the  series.  He  found  that  the 
accuracy  is  a  function  of  n.  Bekker's  solution  yields  the 
following  equation: 
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3  W 

(  3—  n ) ( k  ^  +  bk^ )  v  D 


2  n  +  1 


L 


Since  n  <,  2,  according  to  actual  measurements,  n  =  3  is 
excluded  by  practical  considerations.  Erhlich's  equation 
reads  as  follows:  2 


. w  ] 

2n+i  r(]_4 

51  n+.22n2)  - 

#( .  25 

26n-. I4n2) 

0 1 

bk  {TPj 

L 

u 

(bk 

W 

w; 

51n  +.220^ 

- ) 

-  ] 

-2 

2n+  1 

.  .  .  .  8 

Equation  1 

,  may  be  rewritten 

for 

a  driven 

wheel 

as  follows: 

(Figure  5) 

• 

4 

2ti 

2n 

s 

dT  cos  0 

+ 

(  dN  sin 

0  -  R 

=  0  ...  (9 

-a) 

1 

CM 

0 

2n 

-A 

2tt 

52" 

w  +  5 

dT 

si  n  0  - 

dN  cos  0  ; 

=  0 

.  .  .  .(9r- 

b) 

2tt- 

h 

2n 

-A 

M  - 

§  (2tt  dT  = 

0  .  . 

•  •  •  « 

.  .  ..  (9- 

c) 

*-po 
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As  the  wheel  moves  from  position  1  to  position  2,  while 
turning  by  an  angle  0  point,  "P”  is  transferred  to  "P’",  Fig. 
6.  The  path  of  the  point  is  called  a  cycloid  and  its  para¬ 


metric  equation  is  the  following: 

x  =  (9  -  sin  0) . (10a) 

y  =  ( 1  -  cos  0)  . (10b) 


Equation  10  is  valid  if  no  slip  occurs  between  the  hori¬ 
zontal  plane  (  x  axis)  and  the  wheel  perimeter.  In  case  of 


slip  PP"  is  no  longer  equal  to  (D/2)8. 

The  relative  velocity  between  the  fixed  x  axis  and  the 
bottom  of  t^he  wheel  (vg)  determines  the  slip  as  follows: 


(11) 


It  can  be  seen  that  i  is  positive  if  v  and  v„  are 

o  si 

of  opposite  sense  (driven  wheel).  The  theoretical  velocity 
(vT)  is  defined  by: 


D  _  D  e 

2  u  -  2  t 


(12) 


Thus? 


PP”  = 


£  ©  +  v  t 

c.  S 


D 


2  0  "  io  VT  t 


“  ©  (1  -  i  ) 
2  0 


(13) 
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Thus,  the  equation  of  the  path  becomes! 


x 


y 


Denote 

r 


Then 

x 

y 


|  [0  ( 1  -i  )  -  sin  0]  . (14a) 


|  (1  -  cos  0)  . . ( 14b) 

^  (1  -  1  )  . (15) 

2  o 

rO  -  ^  sin  0  . (16a) 

|  (1  -  cos  0)  .  *  . . ( 1 6b ) 


or  using 

0  =  ui 

x  =  ^  [wt(l  -  i  Q )  -  sin  (ut)] . (17a) 


y  =  D/2  (1  -  cos  ut) . .  •  (17b) 

x  =  rut  -  sin  <wt) . (18a) 

y  =  D/2  (1  -  cos  ut')  . .  ,  .  (10b) 


y  =  D/2  (1  -  cos  ut)  . .  ,  .  (10b) 

The  parametric  equation  of  the  velocity  vector  is: 


vx  =  X 

d  x 
~  dt 

= 

1“ 

[(l  -  i^)  -  cos  Cut)]  . 

.  (19a) 

.  >> 

1! 

> 

_  dX 

“  dt 

D 

2  u 

sin  ut  ....... 

.  (19b) 

Or 

Vx  = 

u[r 

D 

2 

cos 

.(20) 
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When 


G 


ut 


0  or  2  n 


vx  =  u(r  -  D/2)  =  -u  D/2  iQ  =  -vT  iQ  =  vg  ...  (21a) 
vy  =  0  .......  (21b) 

Thus  at  the  bottom  of  the  wheel,  the  velocity  is  equal  to  v  . 

s 

The  speed  or  the  absolute  value  of  the  velocity  is: 

I  V  |  =  J  X2  +  y2  -  (J  [  -1  sin2  t)(l  -  iQ)+  iQ2]  .  .(22) 

For  lo  =  1:  f  v  I  =  |  w 

For  iQ  =  0:  |  v  |  =•■  DU' sin  (j  t.)  . (23) 


Imagine  that  the  entire  xy  plane  is  rigidly  attached  to 
the  wheel.  There  exists  a  point  in  the  plane  where  velocity 
is  zero.  This  point  is  called  the  instantaneous  center  of 
motion  (C).  It  can  be  seen  that  C  must  lie  on  the  vertical 
line  of  symmetry.  (Equation  21).  Every  point  in  the  revolv¬ 
ing  and  translating  plane  may  be  thought  of  as  being  on  the 
perimeter  of  a  wheel  of  radius  (D/2)*.  From  Equation  20: 

=  w'[r  -  (D/2)  J  =  0 


when  (n/2) 


r  -  0/2(1  -  i  ) 
o 


(24) 


Thu  s 


yc  22  D/2  -  0/2(1 


i  )  =  D/2  i 

o  o 


(25) 
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Hence,  r  =  (1  -  iQ)  D/2,  is  the  radius  of  an  imaginary 

wheel  which  is  attached  to  the  actual  wheel  and  has  no  slip 

at  its  bottom,  which  is  at  C.  D / 2 ( 1  -  i  )  =  r,  is  often 

o 

called  the  rolling  radius. 

In  case  of  driven  wheels,  v  <0.  Thus,  i  >0,  r  <D/2, 

S  0 

and  y  >0,  which  means  that  the  instantaneous  center  is  above 
c 

the  x  axis.  For  i  =  1  (100%  slip),  r  =  0,  or  y  =  D/2.  Thus 

o  c 

C  coincides  with  the  center  of  the  wheel.  For  towed  wheels, 

v  >0,  i  <0,  r  >D/2  and  y  <0 
s  o  c 

since  r  =  D/2(i  +  |  iQ|  ^ . (26) 

C  is  below  the  x  axis. 

For  iQ  =  1(100%  negative  slip) 

r  =  ?  D/2  =  D 

Consequently,  iQ  =  -100%  does  not  represent  a  complete¬ 
ly  blocked  wheel.  When  i -«=,  vT  -*■  o,  r  -*■  <»,  or  -  », 

which  means  that  the  wheel  slides  on  the  ground  without  turn¬ 
ing.  (When  vT — ►  a,  u— *o,  since  vT  =  ut).  It  is  not  advis¬ 
able  to  use  two  different  definitions  for  the  slip,  as  is 
generally  done,  d<  pending  on  whether  one  deals  with  driven 
or  towed  wheels,  because  the  mathematics  cannot  be  kept  com¬ 
pletely  general,  In  other  words,  if  one  defines  negative 
slippage  so  that  it  becomes-100%  when  the  wheel  is  blocked, 


so  that 
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then  the  case  of  towed  wheel  < 
of  equations  derived  for  driv< 
An  interesting  geometric 
vector  distribution  along  the 
in  the  following.  Construct  ; 
the  abscissa  of  its  center  be 
its  ordinate  be  zero.  Figure  ' 


The  velocity  vector  associated  with  an  arbitrary  point 


A  on  the  wheel  perimeter  can  be  found  as  follows* 

a.  Connect  A  and  C 

b.  Draw  a  line  perpendicular  to  AC  through  0 

The  line  segment  oW  represents  the  velocity  vector  at 
A. 


Proof t 

-  (ru  +  D/2  w  cos  ^ 

)  i  +  R  u  sin^j 

cos  j 

3  -  -cos  $  -  -cos  0 

sin  ^ 

3  =  sin  0 

Hence  v.  =  (ru  -  D/2  u  cos  9)  T  +  D/2  u  sin  0  J  «  ,  (J 


Equation  27  is  equal  to  equation  19  or  20,  Q.E.D.  (To 
see  that  0  =  consider  the  triangles  AO'C  and  00"A*. 

Two  of  their  sides  are  proportional,  D/2,  r  and  D/2  r  w 
and  one  of  their  angles  is  equal. Hence  the  triangles  are 
similar,  thus  all  three  angles  are  equal). 

It  can,  therefore,  be  concluded  that  the  circle  con¬ 
structed  is  the  hodograph  of  the  path. 

DRIVEN  WHEELS: 

The  wheel  perimeter  may  be  divided  into  four  regions 
according  to  the  "behavior''  of  the  velocity  vectors. 

Region  K  The  velocity  vectors  point  up  and  to  the 
right  along  Section  41  (Fig.  8). 
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Region  21  The  velocity  vectors  point  downward  and  to 


the  right  on  Section  12  Therefore,  if  the  sinkage  of  the 
wheel  zq  is  such  that  zq^  positive  horizontal  deformation 
and  compaction  is  imparted  on  the  soil  along  the  section  as¬ 
sociated  with  zq  -  yc«  The  point  on  the  wheel  perimeter 
whose  ordinate  is  y  belongs  to  an  angle  0  =  cos'd  -  i  ), 

C  VO 


because  (from  Equation  19a)  v 


when  ut  =  0  =  cos'd  -  i  ). 

o 


Thus,  the  region  characterized  by  z  -  y  may  also  be  defined 

0  c 


as  follows! 

2  rr  - 

or  from  Equation  3! 


f- 


£0^0 


2  z 


2  tr  “  cos'd - jp -)  <  0  <,  cos'd  -  io> 

The  tangential  forces  acting  on  the  wheel  in  Region  2 
have  negative  horizontal  components,  thus  an  additional  re¬ 
sistance  will  occur.  The  sum  of  these  tangential  forces  has 
been  denoted  (38).  Note  that  Region  2  is  important  only 


when 


2  IT 


‘A 


„  s  ®»  =  coil(1  -  ‘o’ 


Thus,  the  sinkage  and  the  slip  has  to  satisfy  the  above 
inequality  when  H 2  ^  0. 

Region  3!  The  velocity  vectors  point  downward  and  to 
the  left  along  Section  23  Compaction  and  positive  horizon¬ 
tal  shear  stress  components  occur  here  when  the  wheel  is 
driven.  The  section  of  the  wheel  perimeter  along  which 
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positive  shear  stresses  (tractive  forces)  occur  is  defined 

by  BQ,  When 
• 

arc  23 

Region  Here  the  velocity  vectors  point  upward  and 

to  the  left.  Neglecting  the  small  recovery  effect  of  soils 
this  region  is  of  no  importance,  since  the  wheel  is  not  in 
contact  with  the  ground  along  arc  34 L 

Zero  Slip!  For  zero  slip,  r  =  D/2(l  -  iQ)  =  D/2.  The 
hodograph  is  tangent  to  the  y  axis.  (Fig.  9).  (C  is  at  0). 
There  are  no  velocity  vectors  with  negative  horizontal  com¬ 
ponents.  Hence,  traction  cannot  be  developed. 


Fig.  9 


2  TT  - 


r- 


<,  0  this  section  becomes  the  entire 
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There  are  two  regions  only.  The  velocity  vectors  along  OB 
point  upward  and  to  the  right,  while  along  BO  the  vectors 
point  downward  and  to  the  right. 

Negative  Slip!  For  negative  slip,  the  diagram  is  shown 
in  Fig.  10.  There  are  two  regions  ns  in  the  no-slip  case, 


but  v  =  v  at  the  bottom  of  the  wheel. 


Acceleration! 


The  acceleration  of  a  point  located  on  the  wheel  per¬ 
imeter  may  be  obtained  from  Equations  19a  and  b. 

2 

a  =  x  =  =  dvx  =  D/2  u^  sin  ut  .  .  (28a) 

dt^  dt 

d2v  dvv  ,  2 

a  =  v  =  — 5-  =  -rp-  =  D/2  u  cos  ut  ...  (28b) 

y  j  dt  ax 

Note  that  the  acceleration  is  independent  from  the 

slip.  It  can  be  seen  that  the  endpoints  of  the  acceleration 

2 

vectors  describe  a  circle  of  radius  D/2  w  .  The  magnitude 
of  the  acceleration 

2 

=  D/2u  is  constant. 

Its  direction  is  parallel  to  the  radius  and  its  sense  points 
toward  the  center  of  the  wheel,  just  ,»s  in  the  case  of  a  uni¬ 
form  circular  motion. 


The  radius  of  curvature  is! 

9  ,  3/2 

u  +  r> 

♦  * 

x  y 

*  •»  •  « 

X  y 

For 

iQ  0  ^  =  2  D  sin  (--  t) 


( 29 ) 


(  20  ) 


T  h  u  s 


(>  ~  2  AC 


0 


for  i  = 
0 


LaTl  =  D/2  u.2  cos ( 2"  t) 


.  (33) 


for  1  =  1 


I  aT/ 


=  0 


The  normal  component  of  the  acceleration  is 


J;  |  =  li!i  -  “  u 

•aN 1  o  "2 


2  sin2(£  t)(l  -  i  )  +  i 
2  00 

[4  sin2(^  t ) ( 1  -  i  )  +  i  2] 
2  00 


(34) 


for  i  =  0 
0 


(_  I  D  2  u 

aN I  =  2  “  sin  <2  U 


(35) 


for  i  =  1 
o 


I-  J  _  D  ,2  jvj2 

JaNl  2  u  D TT 


From  Equations  32  and  34 


1*1  =  \N  +  N 


D 


2 


. (36) 


This  checks  with  Equations  28a  and  b.  Thus,  when  the  wheel  moves 
with  constant  speed  u  =  const)  the  acceleration  vectors  are  not 
different  from  that  of  uniform  circular  motion. 

SHEAR  STRESS  STRAIN  RELATIONSHIP 

The  next  task  is  to  establish  a  shear  stress-strain  re¬ 
lationship  which  allows  one  to  describe  an  experimental  soil 


29 


shear  curve  by  means  of  coefficients,  which  depend  on  the  soil 

and  its  state.  (Moisture  content,  density,  load  history,  etc.). 

(39) 

The  soil  shear  strength  was  first  expressed  by  Coulomb  , 
as  followst 

s^  „  =  c  +  p  tan  t . .  (37) 

max. 

The  numerical  values  of  c  and  4  may  be  obtained  by  tri- 
axial  tests  or  by  direct  shear  test.  The  well  known  triaxial 
test  procedure  is  as  follows^®\  A  cylindrical  soil  speci- 
man  is  subjected  to  p^  axial  stress  and  p^  radial  stress. 

Mohr1 s  diagram.  Figure  12,  represents  the  stresses  act¬ 
ing  on  a  plane  which  inclosed  (-|  -  )  angle  with  the  axis 

of  symmetry. 


Figure  12. 
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When  failure  occurs  the  shear  stress  reaches  the  shear 
strength  in  the  failure  plane,  according  to  Coulomb’s  cri¬ 
terion.  Thus,  s  has  to  satisfy  the  equation  of  Mohr's 

max 

circle  and  that  of  Coulomb's  straight  line  (Equation  37). 


Thus, 


Pi  "  P. 


sin  ^  -  c  +  p  tan  (i 


(30) 


(Here  Pj  and  are  the  stresses  applied  at  failure). 


Since 


i-  4 


+ 


f  =  o(  =  JL  +  £ 

A  C  A  T  ^ 


if  for  the  failure  plane, 


Thus,  can  be  evaluated  when  cA  p  is  known. 

From  Equations  38  and  39! 

P  i  —  P  o  n  . 

~  sin  (-r-  +  jO  =  c  +  p  tan  0  .  .  (40) 


P1  ~  P3 


cos  =  ,c  ,*t  Pp  tan  i 


c  =  pi  :  p3 


cos  -  p  tan  6 


(41) 


The  value  of  p,  however,  is  still  unknown. 
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From  Figure  121 


inclination  is  4  and  its  intercept  with  the  s  axis  is  c. 

Since  a  triaxial  test  is  a  slow  procedure,  it  is  not  suit¬ 
able  to  obtain  a  large  number  of  soil  values  in  a  short  time. 
Therefore,  quicker  procedures  have  been  devised.  Bekker  pro¬ 
posed  a  shear  apparatus  in  which  a  normal  load  and  a  torque 
is  applied  on  an  annulus  which  rests  on  the  soil  surface  (17, 
41).  Thus,  the  average  normal  and  shear  stress  acting  in  the 
horizontal  plane  directly  under  the  ring  (sQgv)  can  be  re¬ 
corded.  The  shear  stress  which  is  recorded  at  failure  is 
considered  equal  to  smgx>  =  c  +  p  tan  rf,  where  p  is  Hie  normal 
pressure  applied. 

Actually,  it  is  not  evident  that  the  failure  occurs  in 
the  horizontal  plane.  Thus,  the  recorded  value  may  be  a 
shear  stress  value  which  acts  in  the  horizontal  plane  whije 
failure  occurs  and  smgx.  may  act  in  some  other  plane.  Ex¬ 
perimental  evidence  clearly  shows,  however,  that  the  maximum 
torque  per  annulus  area  reading  is  a  linear  function  of  the 
normal  pressure  applied.  The  relationship  may  be  expressed 
by  means  of  an  equation  similar  to  Coulombfls. 


Figure  14. 


Here  index  "B"  stands  for  "Bekker' s  coefficient" f  as 
opposed  to  Coulomb’ s  cohesion  and  friction  angle.  The 
relationship  between  cQ  and  c,  as  well  as  between  and 
is  not  clear  as  yet.  Since  c^  and  is  more  suitable  to 
describe  vehicle  operation  in  soft  soils  than  c  and  ,  Land 
Locomotion  Mechanics  studies  use  the  former  sets  of  param¬ 
eters.  Thus,  cQ  and  will  be  used  and  index  "B"  will  be 
dropped  henceforth. 

An  experimental  shear  stress-strain  curve  is  of  the 
shape  shown  in  Figure  15  in  most  cases.  Sometimes  a  hump 
and  decay  occurs  as  shown  in  Figure  16: 


Figure  15.  Figure  16. 

Bekker  introduced  an  empirical  equation,  similar  to  that 
derived  for  an  overdamped  one  degree  of  freedom  spring-mass 
system  and  proposed  two  additional  soil  values  (Kj  and  to 
replace  the  damping  and  the  natural  frequency  of  the  undamped 
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(5) 

system  .  Weiss  constructed  a  nomogram  which  enabled  one 
to  evaluate  and  Practice  has  shown,  however, 

that  the  hump  bears  no  practical  importance  because  it  occurs 
only  at  certain  undisturbed  clay  soils  which  are  not  surface 
materi  als  .  Furthermore,  the  effect  of  a  hump  may  be 
neglected  in  the  integration  of  shear  stresses  over  a  solid 
boundary  because  the  deformations  associated  with  the  hump 
occur  only  under  a  relatively  small  part  of  the  soil-vehicle 
interface  surface.  Experiments  show  that  the  hump  decreases 
and  finally  disappears  as  p  increases.  When  the  curve  is 
similar  to  that  shown  in  Fig.  15,  the  following  expression 
lends  itself  to  replace  Bekker* s  equation! 

s  =  s  (1  -  e^K) . (44) 

max 

Equation  44,  cannot  be  derived  by  a  limit  process  from 

Bekker’s  equation  unless  some  approximative  assumptions  are 
(45) 

made.  .  Similar  relationships  were  arrived  at  by  Nuttal 
(13),  Soehne  (29)  and  Weinblum  (26). 

The  physical  meaning  of  K  can  be  seen  in  Fig.  17.  Ac¬ 
cordingly,  K  is  the  abscissa  of  the  intersection  of  the 


curve  (s  =  s  )  and  the  tangent  drawn  at  the  origin 
max 

J  J  -  0 


Since  ^ds ^ 

dJ  j=  o 


max 

K 


the  equation 


(ft) 

dJ  j= 


S 

max 


is  satisfied  when  K  =  j.  Thus,  K  is  the  abscissa  of  the 
point  described  above. 

Therefore,  one  can  obtain  the  numerical  value  of  K  by 
drawing  the  tangent  at  the  origin  and  finding  its  intersect 
with  the  asymptote. 

Reece  suggested  another  method  to  evaluate  the  numer¬ 
ical  value  of  K  in  his  discussion  of  Paper  No.  41  at  the 
First  International  Conference  on  the  Mechanics  of  Soil- 
Vehicle  Systems.  (Also  see  Ref.  45).  He  reasoned  that  if 


s 

Smax. 


then 


log(|  -  1)  = 
max. 


I 

K 


So  the  log  of  — —  -1  is  proportional  to  j.  (The  factor 

smax  ,  s 

of  proportionality  being  K).  Therefore,  if  one  plots| 

I 


in  a  x\ 


\  in 

lyon  a  logarithmic  axis  and  j  on  a  linear  scale,  K  will 

represent  the  slope  of  the  straight  line  j  --  )og(—  -  1) 


m  ax# 
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Equation  44  only  approximates  a  true  shear  curve,  hence 
one  seldom  obtains  a  continuous  straight  line  when  replotting 
the  experimental  shear  curve  on  n  seni-log  paper.  Therefore, 
the  argument  for  Reece's  method  emphasizes  that  one  is  free 
to  consider  a  larger  portion  of  the  curve  than  the  initial 
one  which  is  often  not  clearly  definable. 

Next,  the  question  arises  as  to  how  K  is  influenced  by 
the  normal  load  (p),  the  dimensions  of  the  Bevamcter  annulus 
and  possibly  some  other  factors.  In  other  words,  does  K 
solely  depend  on  the  soil  and  its  state?  It  is  most  likely 
that  this  is  not  so.  A  limited  number  of  test  results  ob¬ 
tained  at  the  Land  Locomotion  Laboratory  seem  to  indicate 
that  K  is  proportional  to  the  circumference  of  the  annulus, 

K  also  increases  with  the  normal  pressure,  but  it  becomes 
constant  when  the  normal  pressure  surpasses  3  -  4  psi. 


J.  Adams^^,  investigated  this  problem  in  his  Master's 
Thesis  completed  under  the  supervision  of  A.  Reece.  He 


found  that  K  is  proportional  to  s 


max. 


following  equation: 


s  =  smov  (1  -  e 
max. 


He  recommended  the 


max. 


Here  K  refers  to  the  tangent  modulus  measuredat  s  =1 

lud X  t 

lb/in  ,  In  his  experiments,  shear  tests  were  carried  out 
by  a  shear  block  and  not  an  annulus  and  the  "semi-log”  tech¬ 
nique  was  used  instead  of  the  "tangent  method"  to  evaluate  K. 

It  is  suggested  that  a  broader  Investigation  Is  required  to 
elucidate  the  true  nature  of  the  various  effects  which  in¬ 


fluence  K. 


INTEGRATION  OF  SHEAR  STRESSES 


The  term 


2tt 

H  =  J  „ 

2r[-Po 


dT  cos  0 


(Equation  9-a) 


is  analyzed  in  the  following.  Clearly: 


dT  =  s  b  |  d  0 


(45) 


I V*' 


Figure  19, 
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or  by  Equation  44 


dT 


s 


max 


(1 


.(46) 


From  Equation  431 

dT  =  (c  +  p  tan  *0(1  -  cj/K)  d0 . (47) 


Thus 

hD  2u  _  i  tv 

H  =  2“  5  (c  +  p  tan  ^)(1  -  eJ cos  0  d0  .  .(48) 


The  next  problem  is  to  relate  j  and  0. 

If  one  approximates  the  wheel  as  shown  in  Fig.  20,  then  the 
following  assumptions  may  be  made. 


The  horizontal  component  of  T  is  a  function  of  the  normal 

pressure  prevailing  at  z  and  the  horizontal  component  of  the 
total  deformation.  The  annulus  sinks  into  the  ground  during 
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a  shear  test  but  only  the  horizontal  deformation  is  recorded. 
The  soil  values  obtained  correspond  to  horizontal  deformation 
components.  Another  question  to  be  cleared  with  respect  to 
this  problem  is  the  following.  The  normal  pressure  is  Vary¬ 
ing  from  zero  to  p.  The  assumption  that  the  shear  stress 
created  is  the  same  as  measured  under  p  =  constant  may  not 
be  true. 

According  to  Figure  21,  the  deformation  at  0  is 


j o  i 


dx= 

2"-  fi. 


d0  .  .  .  .  (49) 

do 


From  equation  14a*. 


Thus, 


=  b  [o 
2 


«  a  fl  C0 

J2  2  J 


-  i  )  -  cos  0]  .  .  .  .  (50) 


)  [(1  -  i  )  -cos  0]  d0 

2 it-  p  0 

r  o 


.  (51) 


When  ^  _i 

2n  -  &  <  Q£2!t~ cos  (1  -  iQ)  -  (Region  2) 


Hence 


J2  =  \  [<*  -  i0H9  -  2tt  +  /?)  -  sin  0  -  sin  /?]  .  .(52) 
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Figure  21 


As  0  becomes  greater  than  cos^ii  -  io)f  Region  3,  the 
direction  of  the  deformation  changes  and  a  negative  j  will 
be  generated  (Figure  22).  In  order  to  maintain  a  negative 
exponent  in  Equation  40,  j  has  to  be  obtained  for: 

2n  -  <2j»-cos*(l  -  i)<  8  <  2n  as  shown  below: 

I  o  —  o  ~ 


2tt-  c  o  s  1  ( 1 


i  > 
o 


-Ztt-cos1(  1 

=  5 

6 


i  ) 

o 


dx  d0 

d0 


(53) 


4  1 


Wheel 


coS^  1— 1 0 ) 


Soil  Level 


C  Vertioal  Tangent 


Acoi1(l-i0) 


Fig.  22 


Jj  =  [(1  -  ioH0v  -  0)  -  sin  0y  +  sin  9]  .  .  .  .  (54) 


WhtiVe 


0  s  2rr  -  cos 1  ( 1  -  i  ) 
v  0 


(55) 


The  sum  of  the  shear  stresses  along  the  soil-wheel  in¬ 


terface  becomes! 


H  =  H.  -  H. 


(56) 
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where 


Hi  = 


bD 

2 


e 


/ 


2"  J  l 

(c  +  p  tan  <0(1  -  e  d  0  . 


(57) 


H2  = 


bJD 

2 


/ 


J  2/K 

(c  +  p  tan  <0(1  -  e  )  d  0  .  .(58) 


2”  -^0 

Equations  56,  57,  58,  54  and  52,  allow  one  to  evaluate 

the  traction  of  a  wheel  when  the  sinkage  is  larger  than  y  . 

c 

In  other  words,  when  the  sinkage  is  deep  enough  to  cause 
soil-wheel  contact  along  that  region  of  the  wheel  perimeter 
(Region  2)  whose  points  move  forward  and  downward.  Using  our 
notations.  Equation  56  is  to  be  used  when  y6>Q^2v.  -  0V. 

When  2it  *  -  0v  only  Region  3  has  to  be  considered. 
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H 

Thus  hn  2n  _  — K 

H  =  — ~  (c+p  tan  #0(1  -  e  )  d9  .  .  .  (60) 

Z'-Po 

Equation  60,  yields  the  tractive  force  when  (See  Equation  3): 

c o s 1  ( 1  -  io)  *  po  =  cos1(l  - 

Before  Equation  56  or  Equation  60  can  be  programmer?  for 
an  electronic  computer,  p  has  to  be  expressed  as  a  function 

of  0. 

Equation  4,  yields  the  expression  desired: 

k  n  n 

p  =  (^—  +  k^)  [(cos  0  -  cos  0q1  tj~]  ,  ,  (4) 

The  foregoing  method  may  be  summed  as  follows: 

1,  Find  z  (Equation  7) 

o 

2,  Find  q  (Equation  3) 

3*  Find  0  =  cos1  ( 1  -  i 

v  0/ 

4.  If  y?o  >  2tt  -  0v  (Use  equations  4,  53,  54, 

57,  58  and  56  to  find  H) 

5.  If  <  2n  -  0 v  (Use  Equations  4,  59  and  60) 

6.  Subtract  R  (Equation  6)  from  H. 


4  5 


The  following  input  data  are  necessary  to  perform  the 
calculations!  kc,  k^,  n  t  c,  0,  K,  D,  bf  W,  iQ.  The 
variable  is  0  naturally. 
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USACDC  Liaison  Office  (SMOTA-LCDC)  .  2 

Sheridan  Project  Office  (AMCPM-SH-D)  .  1 

U.  S.  Naval  Civil  Engineering 
Research  and  Engineering  Laboratory 
Construction  Battalion  Center 

Port  Hueneme,  California  .  1 

Commanding  Officer 
Yuma  Proving  Ground 

Yuma,  Arizona  85364  .  1 

Harry  Diamond  Laboratories 
Technical  Reports  Group 

Washington,  D.  C.  20025  1 
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Defense  Documentation  Center 
Cameron  Station 
ATTN:  TIPCA 

Alexandria,  Virginia  22314  .  10 

Commanding  General 
Aberdeen  Proving  Ground 
ATTN:  Technical  Library 

Maryland  21005  2 

Commanding  General 

Hq,  U.  S.  Army  Materiel  Command 

ATTN:  AMCOR  (TW) . 1 

ATTN:  AMCOR  (TB) . 1 

Department  of  the  Army 
Washington,  D.  C.  20025 

Land  Locomotion  Laboratory  .  2 

Propulsion  Systems  Laboratory  . 5 

Fire  Power  Laboratory . 1 

Track  and  Suspension  Laboratory  .  6 

Commanding  General 

U.  S.  Army  Mobility  Command 

ATTN:  AMCPM-M60 

Warren,  Michigan  48090 . 3 

Commanding  General 
Headquarters  USARAL 
APO  949 
ATTN:  ARAOD 

Seattle,  Washington  .  2 

Commanding  General 
U.  S.  Army  Aviation  School 
Office  of  the  Librarian 
ATTN:  AASPI-L 

Fort  Rucker,  Alabama  .  1 

Plans  Officer  (Psychologist) 

PP&A  Div,  G3 ,  Hqs ,  USACDCEC 

Fort  Ord,  California,  93941 . • . 1 

Commanding  General 

Hq,  U.  S.  Army  Materiel  Command 

Research  Division 

ATTN:  Research  and  Development  Directorate 

Washington,  D.  C.  20025  . 1 
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Commandant 
Ordnance  School 

Aberdeen  Proving  Ground,  Md . .  1 


British  Joint  Service  Mission 
Ministry  of  Supply 
P.  0.  Box  680 
Benjamin  Franklin  Station 
ATTN:  Reports  Officer 


Washington,  D.  C . 2 

Canadian  Army  Staff 
2450  Massachusetts  Avenue 

Washington,  D.  C.  . . 4 

British  Joint  Service  Mission 
Ministry  of  Supply  Staff 
1800  K  Street,  N.  W. 

Washington,  D.  C . ....6 

Director 

Waterways  Experiment  Station 

Vicksburg,  Mississippi . ...3 


Unit  X 

Documents  Expediting  Project 
Library  of  Congress 
Washington,  D.  C. 


Stop  303  . . .4 

Exchange  and  Gift  Division 
Library  of  Congress 

Washington,  D.  C.  20025 . .  .....1 


He adquarters 

Ordnance  Weapons  Command 
Research  &  Development  Division 
Rock  Island,  Illinois 


Attn:  ORDOW-TB . . 2 

Army  Tank-Automotive  Center 

Canadian  Liaison  Office,  SMOTA-LCAN.  ...  .  4 

United  States  Navy 

Industrial  College  of  the  Armed  Forces 
Washington,  D.  C. 

Attn:  Vice  Deputy  Commandant . .10 

Continental  Army  Command 

Fort  Monroe,  Virginia  .......  .  ...1 
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Dept,  of  National  Defense 
Dr.  N.  W.  Morton 
Scientific  Advisor 
Chief  of  General  Staff 
Army  Headquarters 

Ottawa,  Ontario,  Canada  .  1 

Commanding  Officer 
Office  of  Ordnance  Research 
Box  CM,  Duke  Station 

Durham,  North  Carolina  .  3 

Chief 

Office  of  Naval  Research 

Washington,  D.  C . 1 

Superintendent 

U.  S.  Military  Academy 

West  Point,  New  York 

Attn:  Prof,  of  Ordnance . 1 

Superintendent 
U.  S.  Naval  Academy 

Anapolis,  Md . 1 

Chief,  Research  Office 
Mechanical  Engineering  Division 
Quartermaster  Research  &  Engineering  Command 

Natick,  Massachusetts . .  .  .  .  1 

Mr.  M.  Bahn 
Department  6101 
Chrysler  Defense  Engineering 
Box  1316 

Detroit  31,  Michigan  .  1 
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